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Abstract. We solve explicitly the shifted wave equation 

d?u(x,t) = {A x +^-)u(x,t) 

on Damek-Ricci spaces, using Asgeirsson's theorem and the inverse dual Abel trans- 
form. As an application, we investigate Huygens' principle. A similar analysis is carried 
out in the discrete setting of homogeneous trees. 

1. Introduction 

In the book [T7] Helgason uses Asgeirsson's mean value theorem (see Theorem II. 5. 28) 
to solve the wave equation 



(1) 



dfu(x,t) = A x u(x,t), 

u(x,0) = f(x), d t \ t=0 u(x,t) = g(x). 



on Euclidean spaces M. d (see Exercise II.F.l and its solution pp. 574-575) and the shifted 
wave equation 

' dlu{x,t) = (A x + [^f)u{x,t), 
u(x,0) = f(x), d t \ t=0 u(x,t) = g(x), 



(2) 



on real hyperbolic spaces if d (R) (see Exercise II. F. 2 and its solution pp. 575-577). In 
this work we extend this approach both to Damek-Ricci spaces and to homogeneous 
trees. Along the way we clarify the role of the inverse dual Abel transform in solving the 
shifted wave equation. 
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Recall that Damek-Ricci spaces are Riemannian manifolds, which contain all hyper- 
bolic spaces if d (R), H d (C), H d (M), H 2 (0) as a small subclass and share nevertheless 
several features with these spaces. Before jTTj the shifted wave equation on LL d (R) 
was solved explicitly in [2l| Section 7]. Other hyperbolic spaces were dealt with in 
[TUl [T9| 120] and Damek-Ricci spaces in [25] . All these approaches are awkward in our 
opinion. On one hand, [23], [10] and [19j |20] rely on the method of descent i.e. on shift 
operators, which reduce the problem to checking formulae in low dimensions. Moreover 
[TO] involves classical compact dual symmetric spaces and doesn't cover the exceptional 
case. On the other hand, [25J involves complicated computations and follows two dif- 
ferent methods : Helgason's approach for hyperbolic spaces and heat kernel expressions 
P] for general Damek-Ricci spaces. In comparison we believe that our presentation is 
simpler and more conceptual. 

Several other works deal with the shifted wave equation ([2]) without using explicit 
solutions. Let us mention [7] (see also [HI Section V.5]) for Huygens' principle and 
the energy equipartition on Riemannian symmetric spaces of the noncompact type. This 
work was extended to Damek-Ricci spaces in [3] , to Chebli-Trimeche hypergroups in [T3] 
and to the trigonometric Dunkl setting in P, [3] . The nonlinear shifted wave equation 
was studied in [2S1I21E], first on real hyperbolic spaces and next on Damek-Ricci spaces. 
These works involve sharp dispersive and Strichartz estimates for the linear equation. 
Related LPs-LP estimates were obtained in [21] on hyperbolic spaces. 

Our paper is organized as follows. In Section 2, we review Damek-Ricci spaces and 
spherical analysis thereon. We give in particular explicit expressions for the Abel trans- 
form, its dual and the inverse transforms. In Section [3] we extend Asgeirsson's mean 
value theorem to Damek-Ricci spaces, apply it to solutions to the shifted wave equation 
and deduce explicit expressions, using the inverse dual Abel transform. As an applica- 
tion, we investigate Huygens' principle. Section H] deals with the shifted wave equation 
on homogeneous trees, which are discrete analogs of hyperbolic spaces. 

Most of this work was done several years ago. The results on Damek-Ricci spaces were 
cited in |26j and we take this opportunity to thank Frangois Rouviere for mentioning them 
and for encouraging us to publish details. We are also grateful to Nalini Anantharaman 
for pointing out to us the connection between our discrete wave equation (fT6|) on trees 
and recent works [H [9] of Brooks and Lindenstrauss. 



We shall be content with a brief review about Damek-Ricci spaces and we refer to the 
lecture notes [26] for more information. 

Damek-Ricci spaces are solvable Lie groups S — NxA, which are extensions of Heisen- 
berg type groups iV by A = R and which are equipped with a left-invariant Riemannian 
structure. At the Lie algebra level, 



2. Spherical analysis on Damek-Ricci spaces 
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with Lie bracket 

[(X, Y, z), (X', Y', z')\ = (f X'- fx, z Y' — z'Y + [X, X'}, 0) 
and inner product 

<(x, y, z), (x', r, z ')) = (x, x') wm + (y, r v + * 

At the Lie group level, ^ 

S = 

with multiplication 

(X,Y,z) ■ (X\Y',z') = (X+e^X',Y+e z Y'+±e^[X,X'],z + z'). 

So far N could be any simply connected nilpotent Lie group of step < 2 . Heisenberg type 
groups are characterized by conditions involving the Lie bracket and the inner product 
on n, that we shall not need explicitly. In particular Z is the center of N and m is even. 
One denotes by 

n = m + k +1 
the (manifold) dimension of S and by 

Q = f +k 

the so-called homogeneous dimension of N. 

Via the Iwasawa decomposition, all hyperbolic spaces H d ($t), H d (C), H d (M), H 2 (0) 
can be realized as Damek-Ricci spaces, real hyperbolic spaces corresponding to the 
degenerate case where N is abelian. But most Damek-Ricci spaces are not symmetric, 
although harmonic, and thus provide numerous counterexamples to the Lichnerowicz 
conjecture [13]. Despite the lack of symmetry, radial analysis on S is similar to the 
hyperbolic space case and fits into Jacobi function theory [22] . 

In polar coordinates, the Riemannian volume on S may be written as 5(r)drda, where 



const. 



5(r)= 2 m+1 7r?r(f) 1 (sinhf) m (sinh^ fc 
= 2^§r(f)- 1 (cosh|) fc (sinh|) 



r 



n-l 



is the common surface measure of all spheres of radius r in S and da denotes the 
normalized surface measure on the unit sphere in s . We shall not need the full expression 
of the Laplace-Beltrami operator A on S but only its radial part 

radA= (f ) 2 + {^cothf + |tanh|} 9 



Br 

v ' 



S'(r) 
5{r) 

on radial functions and its horocyclic part 

(3) ^f=m 2 f-Qfj 
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on iV-invariant functions i.e. on functions / = f(X,Y,z) depending only on z. The 
Laplacian A commutes both with left translations and with the averaging projector 

fK r ) = W)f dxf(x), 

JS(e,r) 

hence with all spherical means 

fKr)=sb)[ dyf(y). 

JS(x,r) 

Thus 

(4) (A/)J = (radA)/J. 

Finally A has a spectral gap. More precisely its L 2 -spectrum is equal to the half-line 

(—.-¥]■ ' ; 

Radial Fourier analysis on S may be summarized by the following commutative dia- 
gram in the Schwartz space setting pQ : 

<S(R) 

even 

S(Sy ^ <S(M)cven 

A 

Here 

Hf{\)= [ dxip x {x)f(x) 
J s 

denotes the spherical Fourier transform on S, 

Af{z) = e~% z [ dX [ dYf(X,Y,z) 

the Abel transform, 

Ff(\)= [ dze iXz f(z) 



the classical Fourier transform on K and S(S)^ the space of smooth radial functions 
f(x) = f(\x\) on S such that 



su Pr > o (l + r)^e^|(0V(r)|< + 



oo 



for every M, N e N. Recall that the Abel transform and its inverse can be expressed 
explicitly in terms of Weyl fractional transforms, which are defined by 

/+oo 
d(coshrs) (coshr S -coshrr)^- 1 (-^^) M /( S ) 

for r>0 and for /iG<C, MgN such that Re/z> — M. Specifically, 

A = Cl W^ 2 oW k ) 2 and ^^W^oW^, 
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where c\ = 2 3 2 + n t . More precisely, 

k m 

A f( r ) = ~^ (~ d(coshr)) 2 (~d(cosh§)) 2 / ( r ) 

if n is odd i.e. k is even, and 

f + °° k-l m 

^/( r ) = ^7F/ V cosh/- cosh r ( ds) ( d(coshs)) ( d(cosh§)) / ( S ) 

if n is even i.e. is odd. Similarly, the dual Abel transform 

(5) A*f(r)= (f)\r), where f(X,Y,z) = e$'f(z), 

and its inverse can be expressed explicitly in terms of Riemann-Liouville fractional trans- 
forms IZ^, which are defined by 

Kf( r ) = T(pTM) /^(coahra) (cosh rr - cosh rs)^ (j^ffis) 
for r>0 and for /ieC, MgN such that Re/j>-M. 

Theorem 2.1. The dual Abel transform §} is a topological isomorphism between 
C 00 (R) cvcn and C 00 (S^ = C 00 (R) cvcn . Explicitly, 

A*f{r) = f (sinh|)- m (sinhr)-( fc - 1 ) 7^ /2 { (cosh ^ K% 2 [(sinh j)" 1 /] }(r) 

and 

(^)-V(r) = (^ /2 o^ fc/2+1 ){ (sinh, ) m (smh .)*" l /}(r) 



n-1 



2T-T(2) (n-1)! . 

where c 2 = = . More precisely, 

(^)-V(r) = i±( a ^) f (^) 4 -H(^i) m (^r) fc - 1 /(r)} 
z/n is odd i.e. k is even, and 

MT7(r) = (l^^J' ^t^. ( sinh § T ( sinh *)*/(*) 

if n is even i. e. k is odd. 

Proof. Everything follows from the duality formulae 

/ drAf{r)g(r) = / dxf{x) A*g(x) , 
Jm. J s 

r+oo r+oo 

/ d(coah.Tr)W*f{r)g(r) = / d(cosh rr) f(r) U^g^r) , 
Jo Jo 

and from the properties of the Riemann-Liouville transforms, in particular 

Kl /2 : /C°°(R) even r i+1 C°°(R) even 
for every integer £>—!. □ 
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Remark 2.2. In the degenerate case m = 0, we recover the classical expressions for real 
hyperbolic spaces H n (M) : 

n-l f+OO 

Af(r) = ^& / d(cosh s) (cosh s - cosh r) (n " 3)/2 /(s) , 

2 J r 

*A*f(r) = c-s (sinhr) ~( n ~ 2 ) / ds (cosh r — cosh s)~ 5 ~/(s) , 

2 2 ^ i r(f) (n-2)! 
w/iere c 3 = — = -^r 



v^r(V) 2 ^r(=fi)2 



^- 1 / (r) = (2 7 r)-^ i (- - )^/( r ) 



n-l 



d(cosh r) , 

q — .) - in — 1 \ i n — 3 

(^)-V(r) = * (^P" {(sinhr)- 2 /^)} 



i/ n odd and 



r+oo n 
•A f ( r ) = ~n=T~n / Vcoshs-coshr (~ dl) (~ dfcosh s) ) ^ 



r 



z/ n zs even. 

3. Asgeirsson's mean value theorem and the shifted wave equation 

on Damek-Ricci spaces 

Theorem 3.1. Assume that UeC co (S x S) satisfies 

(6) A x U{x,y) = A y U(x,y). 

Then 



(7) / dx' / dy' U(x\ y') = dx' dy' U(x', y') 

JS(x,r) JS(y,s) J S(x,s) J S(y,r) 

for every x,y^S and r,s>0. 

The proof is similar to the real hyperbolic space case [T7J Section II. 5. 6] once one has 
introduced the double spherical means 

U!i(r, s) — tt~t f dx'^-J dy' U(x\ y') 



x,yVi J ) 8{r) / ^ S(s) , 

JS(x,r) JS(y,s) 

and transformed ([6]) into 

(rad A) r E/J;«(r, s) = (rad A) s C/|;» (r, s) . □ 



Asgeirsson's Theorem is the following limit case of Theorem 13. 1[ which is obtained by 
dividing ([7]) by S(s) and by letting s — > 0. 
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Corollary 3.2. Under the same assumptions, 



dx' U(x',y) = / dy' U(x,y') . 

S(x,r) JS(y,r) 

Given a solution u<EC°°(S xR) to the shifted wave equation 

(8) d*u(x,t) = (A x +^)u(x,t) 

on S with initial data u(x, 0) = f(x) and d t \t=ou(x,t) = 0, set 

(9) U (x, y) = e^' t), 

where t is the z coordinate of y. Then ([HD satisfies ([SD, according to (jSJ). By applying 
Corollary 13.21 to ([2]) with y = e and r= we deduce that the dual Abel transform of 
1 1 — y u(x, t), as defined in ([5]), is equal to the spherical mean /J(|£|) of the initial datum /. 
Hence 

u(x,t) = (AT 1 (f!)(t)- 
By integrating with respect to time, we obtain the solutions 

u(x,t) = f ds (AT\st)(°) 
Jo 

to (jHJ) with initial data u(x, 0) = and d t \t=ou(x,t) = g(x). In conclusion, general 
solutions to the shifted wave equation 



d 2 t u{%,t) = (A x +^)u{x,t) 
u(x,0) = f(x), d t \ t=0 u(x,t) = g(x) 



(10) 

on S are given by 

u(x,t) = (AT\fi)(t) + / dsiAT'igDis). 



By using Theorem I2.1[ we deduce the following explicit expressions. 
Theorem 3.3. (a) When n is odd, the solution to (jTOl) is given by 

M (x,t) = c 4 |( I ^^) T (^^)" 1 { i ^ / dyf(y) 



+ C4 (d(cosh±)) 2 (d(cosht)) 2 1 smht / dy g(y) 



S(x,\t\) 

c 

S(x,\t\) 



3 m -+- k -I 7i — 1 

with C4 = 2 2 7r 2 



(b) When n is even, the solution to f llOp zs given fey 



u ( Xj t) = c 9(^\2 ( a\V I d f(y) 

V ' I d\t\ Vd(coshf)/ Va(coshi)/ jB(x\t\) y/cosht- cosh d{y,x) 



k-l 



& 6 WVa(cosh|)V Va(cosht); J B (x\t\) y/ cosh t- cosh d(y,x) 



» C5 = 2 2 1 7r 2 . 



<s 
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Remark 3.4. These formulae extend to the degenerate case m = 0, which corresponds 
to real hyperbolic spaces H n (M>) : 
(a) n odd: 

«(*.*) = ^m(m^)^\Ai [ dyf(y) 



+ c 4o(£ht))~ {^Tt dyg(y)\. 

n+1 n— 1 

2 IT 2 . 

(b) n even: 



re-f-1 n— 1 

with cq = 2 2 vr 2 



'Srajjil) a/ cosht-coshc%,a;) 

+ cr sign(t) ( ^y)^ 1 / dy ^— 



J(a; a/ cosh t- cosh d(y,x) 



With Cj = 2 2 7T 2 . 



As an application, let us investigate the propagation of solutions u to the shifted wave 
equation fllOp with initial data /, g supported in a ball B(xq,R). The following two 
statements are immediate consequences of Theorem 13 .31 Firstly, waves propagate at 
unit speed. 

Corollary 3.5. Under the above assumptions, 

supp u c { (x, t) G S | d(x, x ) < + R } . 

Secondly, Huygens' principle holds in odd dimension, as in the Euclidean setting. This 
phenomenon was already observed in [25]. 

Corollary 3.6. Assume that n is odd. Then, under the above assumptions, 

supp u C { (x, t) G S | \t\ — R < d(x, xo) <\t\ + R} . 

In even dimension, u(x, t) may not vanish when d(x, x ) < \t\ — R, but it tends asymp- 
totically to 0. This phenomenon was observed in several settings, for instance on Eu- 
clidean spaces in [27], on Riemannian symmetric spaces of the noncompact type [7J, on 
Damek-Ricci spaces jl], for Chebli-Trimeche hypergroups |14] . ... Our next result dif- 
fers from [7J HH H] in two ways. On one hand, we use explicit expressions instead of the 
Fourier transform. On the other hand, we aim at energy estimates as in [27], which are 
arguably more appropriate than pointwise estimates. Recall indeed that the total energy 

(11) S(t) = K(t)+V(t) 

is time independent, where 

Kit) = \ I dx \d t u{x,t)\ 2 



WAVE EQUATION IN RANK ONE 



9 



is the kinetic energy and 

V(t) = \\ dx (-A x -^)u(x,t) u(x,t) 
J s 

= \\ dx {\V x u(x,t)\ 2 -^-\u(x,t)\ 2 } 
J s 

the potential energy. By the way, let us mention that the equipartition of (fiT]) into kinetic 
and potential energies was investigated in [7] and in the subsequent works [I5J [5J H] (see 
also [TBI Section V.5.5] and the references cited therein). 

Lemma 3.7. Let u be a solution to (flOl) with smooth initial data f, g supported in a 
ball B(xq,R). Then 

u(x,t), d t u(x,t), V x u(x,t) are 0(e~ (Q/2)l ' 1 ) 
for every 16 S and t GR such that d(x, Xq) < \t\ — R—l . 
Proof. Assume t>0 and consider the second part 



k-l 



(12) HM)=(^f(ir dy ^ 



a (cosh |)/ \d(cosht)J J B (xt) A/cosht- cosh d(y,x) 



of the solution u(x,t) in Theorem 13. 31 b. The case t<0 and the first part are handled 
similarly. As B(xq, R) C B(x, t), we have 

g{y) _ / j n , iM 



dy 7 = / dy 

B( 



j, ^ -y/ cosht-coshd(y,:r) Jb(x ,R) V cosh t - cosh c%,x) 

and thus it remains to apply the differential operator 



m h — 1 

1Jt V a (cosh 4 W VafcoshtW 



. a (cosh |)/ Va(coshi), 

to {cosht — cosh<i(y,x)} 2 . Firstly 

( aTcoshl) ) 2 { c °sht — cosh<i(?/,a;)} 2 = const, {cosh t— cosh a;)} 
and secondly 

(a7c^h|)) T { cosht - coshd^x)}" 1 



2 m (cosh |) 2 { cosh £ — cosh x)j 2 , 



for some constants aj. As 



cosht — cosh d(y, x) = 2 sinh t+d ^' x ) s i n h * x e * 



we conclude that D t {cosht — cosh d(y, a;)} 2 and hence t>(x, t) are 0(e - ^*). The deri- 
vatives d t v(x,t) and Va;f(x,t) are estimated similarly. As far as V x v(x,t) is concerned, 
we use in addition that 

sinh.d(y,x) = 0(e*) and | V x d(y, x)\ < 1. 
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This concludes the proof of Lemma 13.71 □ 

Theorem 3.8. Let u be a solution to (flOj) with initial data f,g£ C£°(S) and let R = R(t) 
be a positive function such that 

R(t) +oo 

/, as t — > ±oo. 

R{t) = o(\t\) 

Then 

dx { \u(x,t)\ 2 + \V x u(x,t)\ 2 + \d t u(x,t)\ 2 } 



l{x,e)<\t\-R(t) 

tend to as t— >±oo. In other words, the energy of u concentrates asymptotically inside 
the spherical shell 

{xeS | \t\-R{t) < d{x,e) < \t\+R{t)}. 
Proof of Theorem \3.8l By combining Lemma 13.71 with the volume estimate 

volB(e, \t\-R(t)) x e Qm - R{t)} as t^ioo, 
we deduce that the three integrals 

dx \u(x, t)\ 2 , 



J 



d(x,e)<\t\-R(t) 

/J™ ir ,,(,.-/ 1 1 2 
:( 



dx \ W x u(x,t)\' 

' d{x,e)<\t\-R(t) 



I 



dx \d t u(x, t)\ 2 

(x,e)<\t\-R(t) 



are 0(e and hence tend to as H ±oo. □ 

4. The shifted wave equation on homogeneous trees 

This section is devoted to a discrete setting, which is similar to the continuous setting 
considered so far. A homogeneous tree T = T q of degree q + 1 > 2 is a connected graph 
with no loops and with the same number q+1 of edges at each vertex. We shall be con- 
tent with a brief review and we refer to the expository paper [12] for more information 
(see also the monographs [T6| [15]). 

For the counting measure, the volume of any sphere S(x, n) in T is given by 



5{n) 




if n = 0, 
if neW. 



Once we have chosen an origin OgT and a geodesic u : Z — >T through 0, let us denote 
by | a; | e N the distance of a vertex iGTto the origin and by h(x) G Z its horocyclic 
height (see Figure [T]). 

The combinatorial Laplacian is defined on Z by 

£ z /(n) = /(n)- ^)+/("-D , 
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Figure 1. Upper half-space picture of T 3 

and similarly on T by 

(13) ^/W^W-^rE^/fe)- 

The L 2 -spectrum of £ T is equal to the interval [1 — 7, 1 + 7], where 



7 = ql/2 + q - 1/2 e(0,i). 



We have 



(14) C T f(n) 



on radial functions and 



f/(0)-/(l) ifn = 



Cf(h) = f(h) - ^ f(h-l) --L.f( h + 1) 

(15) =7f^{p/(ft)} + (l-7)/(ft) 

on horocyclic functions. 

Again, radial Fourier analysis on T may be summarized by the following commutative 
diagram 

C°°(R/rZ) even 



5(T) 



.A 



5(Z), 



Here 



W A ) = E ^a(*)/(*) VA. 
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denotes the spherical Fourier transform on T, 

h(x)=h 

the Abel transform and 

WHE^z^/W VAgK 

a variant of the classical Fourier transform on Z. Moreover r = iS(Z) cven denotes 
the space of even functions on Z such that 

sup neN * n k \f(n)\ < +00 V/cGN, 

and S(Tf the space of radial functions on T such that 

sup„ eN * n fc gf \f(n)\ < +00 V/cGN. 

Consider finally the dual Abel transform 

|rr|=n 

The following expressions are obtained by elementary computations. 
Lemma 4.1. (a) The Abel transform is given by 

Am = q^mi) + *f ^ +k m\+2k) 

= J2 + k Z q ^ +k ^ fm+2k) - f(\ h \+ 2k + 2 )} v h ^ 

and the dual Abel transform by 

A*f(n) = -^q-^f(±n) + q —±q-^ V f(±k) 

k has same parity as n 

if neZ*, resp. A*f(0) =/(0). 

(b) TTie inverse Abel transform is given by 

•A- l f\n) = J2 + k Z q-^- k {f(n + 2k)-f(n + 2k + 2)} 

= q^f(n)-(q-l)J2 + k Z q ^ kf ( n + 2k) VnGN 
and the inverse dual Abel transform by 

(AT 1 f{h) = \q^f{h) + \q--*f{l) 

+ \ E^ ^" 2fc+1 ( M " 2* + 2) - - 2k) } 

2 * J V ' 2 * ^0<fcodd</i JK ' 

if /iGN is odd, respectively 

(AT 1 f(0) = f(0) 
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and 

{A*y l f{h) = U^f(h) + !<r*/(o) 



\ EL ^" 2fc+1 {/(* " 2* + 2) - /(* - 2*)} 

^ /(/») - q-^ /(0) 



'0<fc even</i v ' 

if h&W is even. 



(16) 



We are interested in the following shifted wave equation on T : 

7£^w(x,n) = (££-l+7)u(z,n), 
tt(x,0) =f(x), {u(x,l) - u(x,-l)}/2 = g(x). 



As was pointed out to us by Nalini Anantharaman, this equation occurs in the recent 
works [SI [9] . The unshifted wave equation with discrete time was studied in [TT] and the 
shifted wave equation with continuous time in [23] . 

We will solve ffTB"]) by applying the following discrete version of Asgeirsson's mean value 
theorem and by using the explicit expression of the inverse dual Abel transform. 

Theorem 4.2. Let U be a function on T such that 

(17) C T x U(x,y) = C T y U(x,y) \/x,yeT. 
Then 

y y u(x', y ') = y y mxW) 

t-^ x'&S{x,m) ^>y'£S(y,n) ^ x'eS(x,n) ^— ' y'eS{y,m) 

for every x,y£T and m,nGN. In particular 

(18) y U(x',y) = y U(x,y'). 

V ' ^x'eS(x,n) K ' ^y'eS(y,n) V y ' 

In order to prove Theorem I4.2[ we need the following discrete analog of (j4j) . 
Lemma 4.3. Consider the spherical means 

f|( n )= i V" f(y) VxGT, VnGN. 

Then 

(£ T f)l(n) = (rad£) B /J(n), 
where rad£ denotes the radial part (jl4"|) of C T . 

Proof of Lemma \4-3\ We have 

/(x)-/J(l) ifn = 0, 

/J(n)-5TT/J(n-l)-5ii/J(n + l) if n.eN*. □ 

Proof of Theorem Fix y £ T and consider the double spherical means 



:£ T /)*' 



y ' x'&S(x,m) v ' y'eS(y,n) 
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that we shall denote by V(m, n) for simplicity. According to Lemma fl~3| our assumption 
(fTTj) may be rewritten as 

(19) (rad C) m V(m, n) = (rad C) n V(m, n) . 
Let us prove the symmetry 

(20) V(m, n) = V(n, m) Vm,nGM 

by induction on £ = m + n. First of all, (I2"0"j) is trivial if £ = and (I2(jp with £ = 1 is 
equivalent to f lT9]) with m = n = 0. Assume next that £ > 1 and that fj20|) holds for 
m + n < £. On one hand, let m > n > with m + n = £ + 1 and let 1 < k < m — n. We 
deduce from fTT9|) at the point (m — A;, n + k — 1) that 

Wm — k + l,n + k — 1) — Vim — k, n + k) = 

(21) 

v ' = q{V(m-k,n + k-2)-V(m-k-l,n + k-l)}. 

By adding up f[2T|) over fc, we obtain 

(22) V(m,n) - V{n,m) = q{V(m-l, n-1) - V{n-l, m-1)} , 

which vanishes by induction. On the other hand, we deduce from ( TT9l at the points (£, 0) 
and (0, £) that 

' V(£ + l,0) = (q + l)V(£,l)-qV(£,0) t 
V(0,£ + 1) = (q + 1) V(l,£) - q V(0,£) 

Hence V(£+1,0) = V(0,£+l) by using ( 122^) and by induction. This concludes the proof 
of Theorem 14.21 □ 

Let us now solve explicitly the shifted wave equation (TT6]) on T as we did in Section 
|3]for the shifted wave equation (fTUj) on Damek-Ricci spaces. Consider first a solution u 
to ( Tl6|) with initial data 0) = f(x) and 1) — u(x, — 1)}/2 = 0. On one hand, 

as (x,n) i — ^ m(x, — n) satisfies the same Cauchy problem, we have u(x, —n) = u(x,n) by 
uniqueness. On the other hand, according to (IT5|) . the function 

h(y) 

U(x,y) = q 2 u(x,h[y)) Vi,?/6T 

satisfies ffTTl) . Thus, by applying ffTSl) to [/ with y = 0, we deduce that the dual Abel 
transform of n 4 u(x,n) is equal to the spherical mean f£(n) of the initial datum /. 
Hence 

u(x,n) = (X) _1 (/i)(n) ViGT, Vn6N. 

Consider next a solution w to HTB]) with initial data w(x, 0) =0 and {m(x, 1) — m(x, — 1)}/2 
= g(x). Then u(x,n) is an odd function of n and 

/ \ u(x,n + 1) — u(x,n — 1) 

v(x,n) = - Lj 2 

is a solution to HTB]) with initial data v(x,0) = g(x) and {v(x, 1) — v(x, — 1)}/2 = 0. 
Hence 

2 > v(x,k) if neN* is even, 

-^0<fcodd<n v ' 

5f(x) + 2^ u(x, k) ifnGN*isodd, 

< * < k even < n 



u(x, n) 



WAVE EQUATION IN RANK ONE 



15 



with v(x,n) = (A*) 1 (gl)(n). By using Lemma H~Tl b. we deduce the following explicit 
expressions. 



Theorem 4.4. The solution to ffT6j) is given by 

u{x,n) = §g~»yV . . .f{y) - ^r<r*y2,, ^, , f(v) 

L — < d(y,x) = \n\ z — ' d(y,x)< \n\ 

\n\—d(y,x) even 

+ sign(n) V g(y) VxGT, VnGZ* 

*-~sd(y,x) < \n\ 

\n\—d(y,x) odd 

In ot/ter words, q g 



(23) u(z,n) = M| "' 2 M| "'- 2 /(z) + signHMi^^x), 

n — d(y,x) even 

»/ n>0 and M_i=0. 

Remark 4.5. Notice that the radial convolution operators C n and S n above correspond, 
via the Fourier transform, to the multipliers 

\ t sin nA 

cos nA and . , , 

where cos q X = - — ^ — an( ^ sin g A = 2 — — . 

As we did in Section [31 let us next deduce propagation properties of solutions u to the 
shifted wave equation ffT6l) with initial data /, g supported in a ball B(xq,N). 

Corollary 4.6. Under the above assumptions, 

(a) u{x,n) = 0(q-^) ViGT, VnGZ, 

(b) supp n C { (x, n) G T x Z | d(x, xq) < |n| + iV } . 



Obviously Huygens' principle doesn't hold for (116]) . strictly speaking. Let us show that 
it holds asymptotically, as for even dimensional Damek-Ricci spaces. For this purpose, 
define as follows the kinetic energy 



and the potential energy 

(25) «*(*.w)=a 

= cX (^«-7)«(^n)u(a;,n) 

for solutions u to (TTfi"|) . Here 

#M = /(*) - hffti £„ eSM M 
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is the 2-step Laplacian on T and 

7 = ^(0,1). 

Lemma 4.7. (a) The L 2 -spectrum of C is equal to the interval [7, ]■ Thus the 
potential energy (|25j) is nonnegative. 
(b) The total energy 

£{n) = K(n)+V(n) 

is independent of nGZ. 

Proof, (a) follows for instance from the relation 

C = z±lc T (2-C T ) 

and from the fact that the L 2 -spectrum of C T is equal to the interval [1 — 7, 1 + 7]- 
(b) Notice that the shifted wave equation 

>yClu(x,n) = (C T x -l+j)u(x,n) 

amounts to 

u(x,n+l) + u(x,n-l) = -j= ^2 yeS(x u(y,n). 

As 

E> u(y, n)u(z,n) — (o + l) > \u(x, n)\ 2 + > u(y, n) u(z, n) , 

d{y,z)=2 

we have on one hand 

d(x,y)=2 d(x,y)=l 

On the other hand, 

(27) Hn) = %tf E, eT I n )! 2 -^E , j2/eT U M ^) • 

d{x,y)=2 

By adding up (j2SJ) and (127)1 . we obtain 

*w = 2 E, eT i^ n )i 2 + 1 E, eT i^^ ±:L )i 2 

d(cc,j/)=l 

and we deduce from this expression that 

S(n) = £(n±l). 

This concludes the proof of Lemma 14.71 □ 



WAVE EQUATION IN RANK ONE 



17 



Remark 4.8. Alternatively, Lemma \^7\ b can be proved by expressing the energies JC(n), 
Vin), £{n) in terms of the initial data f, g and by using spectral calculus. Specifically, 

KM = I E ^ C ^- C n-lff(x) W) 

+ I E ^ Sn + 1 ~ s n~i) 2 g(x) g{x) 

+ fReV JCn+i- C„_i)(S n+ i- S n -i)f(x) g(x) 

and 

W = iE PT (i-^)0(^)7R 

+ iE x . eT ( 1 -^)^(x)^y 

+ |Re^ (l-L7 2 )L7 n S n /(x)^). 
i/ere we have used the fact that 

a±i(£-7) = !(3-M 2 ) = J(l-C 3 ). 

i/ence 

£ w = E xeT ^/( x )^) + E xeT ^ + ^)sM + 2 Re E xeT ^ + /( x )^) > 

where 

U£ — I (Cn+1— Cn-l) 2 + f (1— C2) Cj, 

Ki + — I (Sn+i— s n -i) 2 + c 2 ) S^, 

J^n!" = I (Cn+1 — Cn-l) (Sn+l" Sn-l) + j (1— Cq) C n S n . 

By considering the corresponding multipliers, we obtain 

U+=\(l-C 2 ), C=|, Wtf = °. 

and we conclude that 

= I E^ 1 " C 2 )/(x)7R+ I E seT l^)! 2 = £ (°) ■ 

Let us turn to the asymptotic equipartition of the total energy S = S(n). 

Theorem 4.9. Let u be a solution to ffToT) with finitely supported initial data f and 
g. Then the kinetic energy K.(n) and the potential energy V(n) tend both to S/2 as 
n — > ±00 . 

Proof. Let us show that the difference K.(n) — Vin) tends to 0. By resuming the com- 
putations in Remark 14.81 we obtain 

/c(„) - v(n) = j2 xeT u-f(x) tm + E xeT V n9( x ) HG) + 2 Re E xeT w nf( x ) 9® > 
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with 

U n = \ (Cn+1— Cn-l) 2 — \ (1— C 2 ) C% = — \ (1— C 2 ) C 2n , 
Ki = | (Sn+1 - 'S'n-l) 2 — | (I - C2) = I ^ 2n ' 

= I (Cn+1 — C n _i) (5^+1— S'n-l) ~~ 4 (1 ~~ C2) Cn S n = — j (1 — C2) S 2n • 

As 

|| C 2n f\\ t ~<*^q-W H/11,1 and || (1- C 2 )f \\ e x < { q -f^ + 2 } || / , 
the expression 

tends to 0. The expressions 

^2 xeT V n9(x)g(x) = \^^C 2n g(x)g(x) 

and 

Y, x&T W nf^)W)= ~ | E xeT ^n/(^) (1-C a )/(X) 

are handled in the same way. This concludes the proof of Theorem 14.91 □ 
Let us conclude with the asymptotic Huygens principle. 



Theorem 4.10. Let u be a solution to (116]) with finitely supported initial data and let 
(N n ) n£ z be a sequence of positive integers such that 

N n +oo 

,, , N as n-> ±oo. 
x iV n = o(|n|) 

Then the expressions 

ExeT l u ( x ' n )| 2 ' Ex.yeT \ u(x,n)-u(y,n)\ 2 , E« £ t I ^(a:, n + 1) - n(a;, n - 1) | 2 

H<|n|-iV„ |a|,|v|<H-JV„ \x\<\n\-N n 

d(x,y)=2 

tend to as n— )-±oo. In ot/ier words, the energy of u concentrates asymptotically inside 
the spherical shell 

{xeT | |n|-JV n < |x| < |n| + iV n }. 
The proof is similar to the proof of Theorem 13.81 □ 
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